1 



Renormalization of Currents for Massive Fermions* fermilab-conf-98/332-t 

Andreas S. Kronfeld and Shoji Hashimoto a 



iep-lat/9810042 



a Theoretical Physics Group, Fermi National Accelerator Laboratory, Batavia, Illinois, U.S.A. 



00 

On 

> 
o 

> 

^t- 
O 
O 

00 
ON 
-t— > 

Jh 

Oh! 



X 



The renormalization of vector and axial- vector currents for massive fermions (in the "Fermilab formalism" 
discussed. We give results for non-degenerate masses, which are needed for semi-leptonic form factors. 



1. INTRODUCTION 

Like many groups, we and our collaborators are 
calculating form factors of semi-leptonic decays, 
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and similar decays to the vector mesons D*, K* , 
or p. This paper discusses the renormalization — 
including effects of the nonzero quark masses — of 
the weak currents inducing the transitions. The 
form factors are needed to an accuracy of a few 
per cent, so it would be ideal to devise a fully non- 
perturbative renormalization program. We show 
here how certain ratios of transitions induced by 
vector and axial- vector currents take care of most, 
but not all, of the renormalization nonperturba- 
tively. We also give the residual one-loop radia- 
tive corrections to these ratios. They are a few 
per cent. One should expect, therefore, that the 
uncalculated higher-order corrections will not be 
needed for some time. 

Our analysis differs from that of the Alpha 
Collaboration || in its treatment of the quark 
masses. First, in no decay may the quarks be 
considered degenerate. Second, the masses of the 
charmed and bottom quarks are larger than the 
energy scale of chromodynamics, Aqcd- Con- 
sequently, chiral Ward indentities are not help- 
ful for normalizing axial-vector currents. Third, 
these masses are usually, in practice, not small 
compared to the lattice cutoff. Therefore, to iso- 
late the leading mass dependence into coefficient 
functions (rather than matrix elements) we ap- 
ply the mass-dependent improvement program of 
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Ref. H to the clover action. In particular, we 
keep nonzero quark masses in Feynman diagrams; 
for a transition such as b — > u the up quark's mass 
can always be set to zero later. 

At tree level currents suitable for massive 
fermions are given by [Spf 
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Z v ^ c 7 ^ b (3) 
Za^j^o (4) 

where, with n denoting the hopping parameter, 

= V2k7 (l + ad{j ■ ipt ', 
= s/2k7 (V - ad{{D&) ■ 7) , 



(5) 



and ip is the field of flavor / in the hopping- 
parameter form of the action. In the following, 
we speak of charmed and bottom quarks, but c 
and b really stand for distinct quark flavors. 

The renormalization factors Z depend on both 
masses in the current. At tree level the mass de- 
pendence factorizes 
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where M( is the rest mass of flavor /. Beyond 
tree level Zy ^ Za and the mass dependence no 
longer factorizes. The main aim of this paper is 
to obtain the full mass dependence of Zy and Za 
in one-loop perturbation theory. Note that for 
degenerate quarks the vector current's absolute 
normalization can be used to define Zyss non- 
pert urbatively. 

The coefficient d{ depends on the mass of fla- 
vor /. Beyond tree level there is no reason to 
expect that a universal rotation a la Eq. (|^) im- 
proves all currents. One strategy is to define d{ 
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by requiring that the equal-mass vector current 
be conserved. For other currents and when the 
masses are not the same, further improvement 
is attained by adding higher-dimension terms — 
dj^iuij^, VD^, and 9i*7 5 *— to the right- 
hand sides of Eqs. (||) and (||), as is done with very 
light quarks Q] and in nonrelativistic QCD |?],||]. 



2. VECTOR CURRENT 

Most of the renormalization can be captured 
by writing 
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The equal-mass factors Zv« are obtained non- 
perturbatively, leaving the expansion 
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The one-loop term Ry\ b is set by requiring 

to one loop, where 

_ 3Cf / m b + m e m b 
167T 2 I mj, — m c m c 



(9) 
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The only Feynman diagram needed to com- 
pute Ry\ b is the vertex diagram. Everything 
associated with the legs — self energies and the 
factors e Ml / 2 e Ml / 2 — drop out by construction. 
From the vertex diagram with unequal masses one 
must subtract the average of the diagram with 
equal masses. From this combination one also 
must subtract the corresponding combination of 
the continuum vertex diagram. This gives —Ry\ b . 
The resulting loop integral is ultraviolet and in- 
frared convergent. (Infrared cancellation occurs 
point-by-point if the masses in continuum prop- 
agators are taken equal to the corresponding ki- 
netic masses.) 

Before showing results it is worthwhile to an- 
ticipate the outcome. See Fig. |l|. Rycb is sym- 
metric under interchange of b and c. When the 
two masses are equal R y \ b = 0, I > 1, by con- 
struction. When both masses are far below the 



Figure 1. Limiting behaviors of Ry cb , R^ cb - 



lattice cutoff, Ry\ b must vanish as m& iC a to a 
power. (The linear term should be proportional 
to 1— cswj an d it is.) When both masses are much 
larger than the cutoff, the heavy-quark flavor 
symmetry (of the lattice action) ensures that the 



lattice contribution vanishes, i.e., R 
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When one mass is larger than the cutoff and the 
other smaller, R v \ b should resemble the radiative 
correction of a "static-light" current: dependence 
on the light mass should drop out, leaving a log- 
arithm (3Cf/167t 2 ) In M\a plus a constant |J. 

Fig. ^ bears out these expectations for a set 
of masses with m c = 0.256m;,, appropriate to 
charmed and bottom quarks. These results have 
been obtained independently by the two authors. 

To illustrate how to apply these results, let us 
consider the quark masses relevant to our cal- 
culation of the semileptonic decay B — > Dlv at 
/3 = 5.7 (and csw — 1)> for which m\,a — 3.9 

and m c a — 1.0 M. Reading off Fig. one finds 

[11 

Rycb — 0.0076. For the strong coupling we con- 
sider the range a s (7r/a) = 0.19 to a s (l/a) — 0.33. 
Thus, we multiply the bare matrix element with 
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= 1 + 0.0076 x 4tt x a s = 1.025(6). (11) 
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M 2c /M 2b = 0.256 
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Figure 2. Radiative corrections to the b — > c vec- 
tor current. Circles (squares) denote csw = 1 (0). 



Figure 3. Radiative corrections to the b 
axial-vector current. Symbols as in Fig. ^. 



This estimate will be refined when we complete 
the calculation of the BLM 10 matching scale q*. 



3. AXIAL- VECTOR CURRENT 

For the axial- vector current one could compute 
Z A ab / Zyab , with the denominator from Eq. (^) . 
For our analysis of B — > D*lv, however, we need 
the ratio 
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The one-loop term ivfL is set by analogy with 
Eq. (H), replacing Vo and 70 with Ai and 7175, but 



remains as in Eq. 



of R^ Ac b are as for Ry\ 



0|). The limiting behaviors 
Our results are in Fig. || 



4. OTHER COMPONENTS 

Scattering matrix elements of Vi and Aq vanish 
as the three-momenta go to zero. They can take 
the same renormalization factors as Vi and Aq, 
but commensurate effects are the one-loop correc- 
tions to d\ and to the coefficients of other higher- 
dimension improvements. The nonperturbative 
construction of a normalized, conserved VU im- 
plies that perturbation theory is needed only for 
"ra fc - m c " and ll V - A" . 



5. SUMMARY 

By computing most of the renormalization non- 
perturbatively, it is possible to reduce the one- 
loop corrections. This is achieved by requiring 
a correctly normalized, conserved vector current, 
or by extracting physics from a ratio of corre- 
lators, such that Eq. ( |i~2|) applies. Indeed, the 
residual coefficients shown in Figs. || and || are 
small. Hence, the uncalculated higher-order ra- 
diative corrections to semi-leptonic form factors 
are unlikely to be large. 

Fermilab is operated by Universities Research 
Association Inc., under contract with the U.S. 
Dept. of Energy. 
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